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Abstract 

Recent evidence for the top mass in the region of 160 GeV for the first time provides 
an opportunity to use the full power of relativistic quantum field theoretical methods, 
available also for weakly bound systems. Because of the large decay width F of the top 
quark individual energy- levels in "toponium" will be unobservable. However, the potential 
for the tt system, based on a systematic expansion in powers of the strong coupling 
constant can be rigorously derived from QCD and plays a central role in the threshold 
region. It is essential that the neglect of nonperturbative (confining) effects is fully justified 
here for the first time to a large accuracy, also just because of the large F. The different 
contributions to that potential are computed from real level corrections near the bound 
state poles of the tt-system which for F 7^ move into the unphysical sheet of the complex 
energy plane. Thus, in order to obtain the different contributions to that potential we 
may use the level corrections at that (complex) pole. Within the relevant level shifts we 
especially emphasize the corrections of order O(afmj) and numerically comparable ones 
to that order also from electroweak interactions which may become important as well. 
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1 Introduction 



Perturbative expansions in the coupling constant in quantum field theory possess two 
types of applications, the calculation of scattering processes and the computation of pro- 
cesses involving weakly bound systems. Many of the successes of quantum electrodynamics 
(QED) are, in fact, related to positronium, i.e. to the second one of the aforementioned 
applications. The proper starting point for any bound-state calculation in quantum field 
theory is an integral equation, comprising an infinite sum of Feynman graphs. The Bethe- 
Salpeter (BS) equation fulfills this task, and it is well known that in the limit of binding 
energies of O(a^m) the Schrodinger equation with static Coulomb attraction is obtained. 
The computation of higher order corrections to the Bohr levels - beyond 0{a^m) — , how- 
ever, turned out to be far from trivial. It was recognized, though, relatively late that, 
at least conceptually, substantial progress with respect to a systematic treatment results 
from a consistent use of a perturbation theory geared to the original BS equation 0|. 
In that manner, at the same time, nonrelativistic expansions as implied by Hamiltonian 
approaches with successive Fouldy-Wouthuysen transformations are avoided. 

Within the BS-technique it is desirable to have an exactly solvable zero order equation 
different from the Schrodinger equation, otherwise e.g. the approximation procedure for 
the wave function lacks sufficient transparency, especially in higher orders. Considerable 
freedom exists for selecting such an equation. Of course, certain corrections included 
already at the zero level are to be properly subtracted out in higher orders. An especially 
useful zero order equation has been discovered some time ago by Barbieri and Remiddi 
(BR equation 0), but also other equations have been proposed 0. Still, one of the most 
annoying features of all bound-state calculations remains the pivotal rule played by the 
Coulomb gauge. In other gauges, e.g. already the (in QED vanishing) corrections 0{a^m) 
of the Bohr levels imply the inclusion of an infinite set of Feynman graphs P], already 
at that order. Only in very special cases, when certain subsets of graphs can be shown 
to represent together a gauge-independent correction, another more suitable (covariant) 
gauge may be chosen. 

In contrast to positronium in QED, the vast literature on bound state problems in 
quantum chromodynamics (QCD) adheres to a description of the quark-antiquark system 
by the Schrodinger equation with corrections 'motivated' by QCD 0. As long as a rela- 
tively small number of parameters suffices for an adequate phenomenological description 
of observed quantum levels, this approach doubtless has an ample practical justification. 
However, again and again certain deviations from such phenomenological descriptions are 
reported §]. Thus also for this reason a return to more rigorous QCD arguments remains 
as desirable as ever - provided Nature offers a "window" where such methods are applica- 
ble. The standard literature on this subject almost exclusively is based on nonrelativistic 
expansions |0 or on the calculation of purely static Coulomb forces 0. Also very often 
potentials with higher order corrections as determined from on-shell quarkonia scattering 
are used [jlO|. A comprehensive account of this usual approach is given in ll^]. However, 
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one cannot hope to push some of the calculations further without a solid theoretical ba- 
sis which is provided by the Bethe-Salpeter formalism alone. Terms not expressible from 
the outset as "potentials", or relevant off-shell effects, which are typical for higher order 
corrections, may be even lost |T^. The importance of such a perturbatively consistent 
approach is supported by the empirical fact that at least in one case, namely the decay 
of S-wave quarkonium, the result of a full BS-perturbation calculation [|13], including the 
QCD corrections to the bound state wave function, yields a result very different from the 



one which considered only the corrections to the quark antiquark annihilation alone |T4 . 

In this connection the relatively large running coupling constant even at high energies 
represents a well known problem, together with the frequent encounter of large coefficients 
in a perturbative expansion. For this reason e.g. problems arise in the comparison of 
the coupling constant as determined from scattering experiments within the minimal 
subtraction scheme (MS), with the coupling constant to be used in a consistent weak 
bound-state approach. The philosophy within our present work will be that the orders of 
magnitude, as determined from ajjg, will be used for estimates but that we shall imply 
a determination of as by some physical observable of the quarkonium system itself (e.g. 
energy levels, cf. the remarks after eq. ([4.3| ) ). In that way delicate correlations of 'genuine' 
orders of as from basically different types of experiments are avoided. 

Of course, the quarkonium system also differs profoundly from positronium because 
of the confinement of quarks and gluons. However, the phenomenological success of the 
nonrelativistic quarkonium model can be explained by the fact that the bound states 
of heavy quarkonia are sufficiently deep in the Coulomb funnel and thus sufficiently far 
away from the confinement part of the potential. Estimates in the early 80-s |T^ of 



nonperturbative effects, describing the 'tail' of confinement by a gluon condensate [|T| 
suggested that only with quark masses well above about 50GeV the importance of such 
nonperturbative effects for low Bohr quantum numbers may decrease sufficiently to make 
perturbative 'field theoretical' level corrections competitive and observable. An improved 
calculation of that correction based already on the BS approach was done some time ago 



17 1 and shows a result which is smaller by about 20% as compared to the estimate based 



upon the Schrodinger equation. Now from CDF |T^ and from high precision electro-weak 



experiments of the LEP collaborations |jT9[ , the mass range of the top quark seems to be 



established to lie in the range ISO-lSOGel^. Thus for the first time a nonabelian bound 
state quarkonium system seems to fulfill the high mass requirement for a genuine field 
theoretical approach. On the other hand, a special feature of a superheavy top is that 
the weak decay t b + W broadens the energy levels for increasing mass mt so that 
above rut ~ lAOGeV individual levels effectively disappear. This, however, by no means 
invalidates the use of perturbation theory for weakly bound systems. E.g in a future e"*'e~ 
collider experiment at the ti threshold a smooth curve - instead of sharp energy levels 
immediately below the threshold - will the object of experimental studies. Nevertheless 
also this curve can be predicted theoretically. Moreover, the large width Tt^b+w >> ^qcd 
even turns into a virtue since it makes confining effects unimportant: the top has no time 
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to hadronize and thus becomes to a rather high accuracy a fully perturbative object, even 
for QCD pTI, 



In contrast to low mass quarkonia this also offers the possibility to extend bound- 
state perturbation arguments to (nonrelativistic) energies E > above the tt threshold 
E = 0. As shown first in the seminal work of Fadin and Khoze |2T| it is possible to in- 
clude the width in the zero order nonrelativistic Greenfunction by the simple replacement 
E E + iV. The Greenfunction G, whose absorptive part expresses the cross section 
for tt production, beside E + iT only depends on the potential. Previous work along this 
line |2^, ^ was based on QCD-suggested and/or phenomenological potentials. The 
aim of our present study is to find a way to replace those potentials by one, derived from 
systematic QCD perturbation theory only. Any ambiguity in our approach is eliminated 
by the requirement that this potential should be identified order by order with the cor- 
responding one for "nonabelian positronium" , evaluated from (real) energy shifts at the 
position of the bound-state singulatity which has moved into the unphysical sheet of the 
complex energy plane. It should be noted that well known arguments [0] for such an 



aproach even guarantee the gauge-independence of our results, despite the fact that, for 
practical reasons, the calculations are performed in the Coulomb gauge. Proofs of gauge- 
independence are valid order by order in perturbation theory. Also for that reason we are 
careful to keep track of these order, meticulously avoiding e.g. the mixing of orders in a 
running coupling constant. 

To 0{mal) for the energy shifts one encounters the well-known vacuum polarization 
effects, present in a system of two fermions with equal masses m amd gauge coupling 
g = y/Airas, if other fields with lower mass ( here the lighter flavours and the gluons ) are 
present. 

At 0{mag) the 'relativistic' corrections, corresponding to the same type of graphs 
as in the abelian case are present here as well. However, even considering only vertex 



corrections alone, a difference to the abelian case was discovered a long time ago |25 . 
In addition typical other nonabelian contributions may appear at this level as well. In 
contrast to scattering processes, the determination of the order to which a certain graph 
contributes in relativistic bound state perturbation theory in each case requires a special 
analysis. Other terms (from QED and weak interaction), incidentally, may be of the same 
numerical order 0{ma^). The Coulomb gauge also entails peculiar additional nonlocal 
interaction terms in the effective action, appearing e.g. in the path integral formulation 
We also investigate the effect of those terms here. 



In section 2 we recall some basic facts about BS-perturbation theory. The BR equation 
for nonabelian weakly bound onium-systems is found to allow a simple generalization for 
wave functions at complex poles. In fact the advantages of the BR equation really become 
more evident above O(a^). However we want to introduce this technique here, because 
it has shown to yield already at least in one case also a very convenient regrouping 
of corrections: In a computation of the bound state effects for the decay width of a 
"toponium" state a Ward -like cancellation of graphs became evident most clearly just in 
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this approach . 

Tree graphs leading among others to the well known 0{mag) corrections are discussed 
in sect. 3. As indicated already above, the one loop vacuum polarization (sect. 4.1) provides 
a correction term 0{ma^) in the nonabelian case. Here we take the opportunity to point 
out that effects from some of the lighter quarks (especially bottom) in the toponium 
system must be treated more carefully than it is usually done by including the quarks 
only in the number of flavours of a running coupling constant. 

In sect. 4.2 we revisit the old calculation of Duncan for the one-loop vertex cor- 
rection within our present formalism, avoiding some approximations made in that early 
computation which allows us even some statements on the term O(a^m). 

Sect. 5 is devoted to an exploration of possible further corrections of the same numer- 
ical order of magnitude as 0{a'^m). We list several candidates of relevant QCD graphs 
and we show in cases of the most simple two loop graphs that such corrections may well 
be relevant. Beside these graphs from QCD, also corrections from the weak interaction 
and QED may turn out to contribute to this order, but the nonlocal Schwinger- Christ-Lee 
type graphs, peculiar to the Coulomb gauge, are irrelevant to 0{a^m) as shown by an 
explicit calculation. Among others we observe that the effect of the Z -Boson has been 
underestimated so far since it gives rise to a singlet triplet splitting which is equally 
important as the usual Breit interaction. 

It should be stressed that our results for the energy shifts - below threshold are - 
of course - in principle applicable to the bottonium system as well. In fact a "rigorous" 



derivation of the corrections in that case has been presented in the last ref. [|10l, taking 
into account perturbatively the gluon condensate as a description of the confinement 
effects. In view of the large magnitude of that correction in the bb system and of the also 
relatively large theoretical error of the latter we feel safer in the tt system where those 
nonperturbative contributions are completely negligible. 

In the last section 7 we summarize our potential for tt in the threshold region, including 
all terms which lead to numerical orders O(ma^) in the level shifts, as discussed in our 
present paper. Still missing pieces to that order are emphasized. 



2 BS-Perturbation of the BR Equation 

A correct formulation of QCD in Coulomb gauge entails not only Faddeev-Popov-ghost 
terms but also the inclusion of nonlocal interaction terms ||2^. Therefore, the full La- 
grangian reads (a=l,...,8 for SU(3)): 

1 ^ 

where the Lagrange multiplier B'^ guarantees the Coulomb gauge, and where 

D, = d,~tgT'^A;, (2.2) 
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F;, = d.At-d^AI^ + gfabcAlAt. (2.3) 



f 1 and f 2 are given in [£6[ and are discussed more explicitely below. The above Lagrangian 
will include all effects of the strong interaction, but, as we will show, QED and weak 
corrections may also give a contribution within the numerical order of our main interest 

mat)). 

The BS equation in terms of Feynman amplitudes for K and S reads for a bound state 
wave function x 

X?^{p;P) = -^Su'{^+p)S,,{~+p) I 0^,K,,,,,,AP,P,P')X^'W-^P)^ (2-4) 

where S is the exact fermion propagator, and K is the sum of all two fermion irreducible 
graphs. Furthermore, we have introduced relative momenta p and p', a total momentum 
P = Pi — p2, and we choose a frame where P = [Pq, 0) = (2m + E,0). The notation can 
be read off the pictorial representation in fig.l. i,j are collective indices for spin (a, p) 
and colour (noted a,/3). 

2.1 Zero order equation 

It is well known that the dominant part in K for weak binding (a ^ 0) is the one- Coulomb 
gluon exchange which results in an ordinary Schrodinger equation with static Coulomb 
potential. This result is even independent of the chosen gauge in the ladder approximation 
by a simple scaling argument po ~ O(ma^), \p\ ~ 0{ma), Pq ^ 2m — 0(ma;^) [p^] : 

K^K,j"e,j°:^-j^^jl.j% (2.5) 

In equation ( p.5|) we have already used the fact that only colour singlet states can form 
bound states because the Coulomb potential is repulsive for colour octets. The colour 
trace will always be understood to be already done, leading to the definition 

6 An 3 

to be used in the following, in terms of the usual strong coupling constant as- Because the 
above mentioned nonrelativistic limit of the BR equation contains the projection operators 
A^, defined below, it is awkward to calculate the so-called relativistic corrections in a 
straightforward way within the framework of BS perturbation theory, starting from (|2.4|) 
with ( |2.5| ). Therefore, we use the BR equation instead of the Schrodinger equation. 

Moreover, for our present case, we need a generalization of 0] for unstable fermions, 
described by complex m — m = m — 'iro/2 where ro/2 represents a suitaby corrected 
(c.f. below) weak decay rate of the free top quark. Since the real part of m should still 
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determine the sum of polarizations in the numerator of the propagator, we take by analogy 
to the stable case in the zero order approximation to ( p.4|) the relativistic free propagator 
{Ep = y/p'^ + m^, T] = mTo/2Ep) 



Si±Y+Po) - [(±Y+P)7-m]-^= (2.7) 

^^70 , A"7o I n(^^ (9^^ 



±1^ +Po- Ep + iT] ±^ + po - Ep - ir] m 
with the relativistic projectors 

AHp) - ^^^(^^"'+^"^). (2.9) 
2Ep ^ ' 

If furthermore dKo/dpo = both sides of ( p.4|) may be integrated with respect to po- On 
the r.h.s. the product of the two propagators S with ( p.7|) yield four terms, with Cauchy 
poles determined by irj, two of which give no contribution. Generically we obtain 

^Pog^g ^ A^7o ® A-70 A-70 (g) A+7o 
2m Pq - 2Ep + 2ir] -Pq + 2Ep - 2ir] ^ ' ^ 

so that Ko{p,p')^o with 

^o = ^[xodpo (2.11) 
remains to be inserted at the places of " (8)" . As in the kernel is now chosen as 

Ko^o = [7oA+A+A+]$o[A-A-A-7o]if (2.12) 
so as to annihilate the second term in ( p.lO|) . In 

{Po-2Ep + 2tri)^o = I ^i^(p,p',Po,^^)A+A+A+<l>oA-A-A- (2.13) 

instead of A~ the projector A^ := A^{—p) appears. Thus expanding $0 = Z1a,_b=± A'^'Io'^'^^A^ 
only "I'o^"'' is found to differ from zero and obeys 

d^p' 
(2^ 

The nonrelativistic projection of (|2.14|) onto eg) with 

A+A+A+ = /iA+ 

A^A^A" = (2.15) 
Ep^ m 



(P„ - 2E, + 2r,)A+4.},+-'A- = ^ i^A'Cp.p', P„,r,)A+A+4.<+-'(p')A-A-. (2.14) 



and the introduction of appropriate factors in K relatively to the Coulomb kernel (|2.5|) 

k = {jjfi'uu'y^K^ (2.16) 
= {Po + 2Ep + 2ir])/4 

and in the wave function 

<l>o(p) oc z//i0(p-) (2.17) 

also for the case of an unstable fermion lead to a Schrodinger equation for the wave- 
functions $(p) in momentum space {E = Pq — 2m + 2iT] = E + 2ir]): 

i^-Km = (2.18) 

The eigenvalues for En = E + ^ in ( p.l^ ) clearly occur at the real Bohr levels En, i.e. 

£„ = M»-2m=-!^-!;^ + 0(a«) (2.19) 

In addition to the selection of the "large" components by the choice ( p.l2| ), obviously 
also the sign of ir] in ( p.l6|) was crucial for the dependence of (|2.18| ) on the combination 
E + 2iT] alone. Still, for the bound-state argument at complex values of the energy, the 
independence of t] with respect to the momentum p is essential. Going back to (|2.7| ) we 
observe that the choice 

f = (2.20) 

m 

with Tt = const the (cm.) decay rate of a single top quark, the disturbing p-dependence 
in 1] is cancelled. We note that ( |2.20|) may be interpreted as the width corrected by time 
dilatation. 

The full BS-wave function (colour singlet, BS-normalized ||2^) for the BR-kernel can 
be obtained by going backwards to Xo^~^ '■= x{p)- It appears as the real residue of the 
complex pole. The wave function (a;„ = Ep — {2m + -E„)/2 — ie) 

I \ -*i \ . A+5Ar /i(p) 2^„ 
Xn(p, e) = 70X„(P, -e)7o = ^ (^gT^^T^'^^^^- (^.21) 

is identical to the BR wave-function of a stable quark and belongs to the spectrum of 
bound states P„ = M° — iVt. The ie has been introduced to determine the integration 
around that pole. 

In eqs. ( [^.21D S" is a constant 4x4 matrix which represents the spin state of the 
particle-antiparticle system: 

S=[l'K, ■ f^f,' (2.22) 
I am") A : triplet. 
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is simply the normalized solution of the Schrodinger equation in momentum space, de- 
pending on the usual quantum numbers (n, /, m) ||30|, a±i, Oq in ( p.22| ) describe the triplet 



states. In the following it will often be sufficient to use the nonrelativistic approximations 
of eqs. (|2.21|) {un ^ oJn = {p^/m + E„)/2 - ie) 

= ::^^^^^S = -^,x{pr^, (2.23) 



The net result of this subsection is that not only in the nonrelativistic case but 
also for the BR-equation, by a suitable choice of the kernel and of the wave function to 
zero order the exact solution can be reduced to the one of the Schrodinger equation at the 
pole with energy E shifted to real values oi E = E + 2iri = E + iTt, where Tt is the decay 
rate of a single top quark in its rest system. Thus the neccesary requirement is fulfilled to 
calculate a rigorous QCD potential from its effects at such a pole. That this pole lies off 
the real axis is unimportant. Another way to express the same fact is the following: Since 
the only effect of the finite width Tf could be absorbed in a (complex) energy alone, also 
for the BR-equation we could perform all calulations at real energies E (or E = E + ie if a 
singularity occurs on the real axis), and continue analytically afterwards to E = E — iTt- 
Clearly in this way the corresponding function in the upper E plane is obtained, because 
the Breit-Wigner pole (in our case ), of course, has a complex conjugate counterpart at 
ImE < 0, seperated by a cut along the real E axis. 

2.2 Level-Shifts 

Having established that - as far as the systematic determination of the QCD potential 
is concerned - we may now just consider perturbation theory for "nonabelian positron- 
ium" , perturbation theory for the BS equation starts from the BR equation for the Green 
function Gbr = Gq of the scattering of two fermions PT 



tGo = -Do + DoKoGo, (2.24) 

which is exactly solvabel. Dq is the product of two zero order propagators, Kq the corre- 
sponding kernel. The exact Green function may be represented as 

-I oo 

G = E X^f-fs irXni' + Greg = Go E (HGoY , (2.25) 

where the corrections are contained in the insertions H. Bound state poles M„ contribute, 
of course, only for Pq < 2m. It is easy to show that H can be expressed by the full kernel 
K and the full propagators D: 

H = -K + Ko + iD-^ -iD^\ (2.26) 
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Since the corrections to the external propagators contribute only to O(a^) ||3^, the per- 
turbation kernel is essentially the negative difference of the exact BS-kernel and of the 
zero order approximation. 

Expanding both sides of equation (p.25|) in powers of Pq — M°, the mass shift is 
obtained 0: 

AM = {ho){l + {h)) + {hogih) + 0{h^). (2.27) 
Here the BS-expectation values are defined as e.g. 

(2^ J ^2.28) 

We emphasize the four-dimensional p-integrations which correspond to the generic case, 
rather than the usual three dimensional ones in a completely nonrelativistic expansion. 
Of course, ( |2.28| ) reduces to an ordinary "expectation value" involving d^p and $(p), 
whenever h does not depend on po and p'q. 

In ( |2.27| ) hi and Qi represent the expansion coefficients of H and Gq, respectively, i.e. 

oo 

H = J2hniPo-My (2.29) 

n=0 

oo 

Go = Y.9n{Po-Mlr-^ (2.30) 

n=0 

Similar corrections for the wave functions are irrelevant in our present work. 



3 QCD Tree Diagrams 

The contributions from the tree diagrams 2. a to 2.c are well-known from the abelian case. 
Fig. 2. a is peculiar for the use of a different zero order equation than the Schrodinger 
equation. It contains the difference between the exact one Coulomb-gluon exchange and 
the BR-Kernel ( ( |2.12| ) at Pn — F^). The exchange of one transverse gluon is represented 
by graph 2.b. The annihilation graph fig. 2.c with one gluon does not contribute in our 
nonabelian case. 

For the sake of completeness and in order to illustrate the present formalism, we exhibit 
first the results for the tree graphs as well. The perturbation kernel for the Coulomb-gluon 
exchange 



is needed for the calculation of the energy shift induced by fig. 2. a using Eqs. ( |2.27|) 
(^), ( p:^ and (^T^ to (PTT^ ). For the spin-singlet we have: 



AM2.a = {H, + Kbr) 
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16^2' ' Pom J (27r)3 7 (27r)3^'^^ ^ (p-p' 



The contribution from the transverse gluon (fig. 2.b) 



- zH^., = ^^7,a^Ln',>A^9^k + (3.3) 

with 

q^p -p, (3.4) 

gives rise to a spin singlet-triplet (magnetic hyperfine) splitting. Because of the 7-' ma- 
trices, the projectors from both wave functions annihilate (p.3|) . This means that two 
factors p7, contained in A^, are needed for a nonzero result which in turn gives rise to 
two extra orders of a. By this mechanism we arrive at the well known contribution O(a^) 
from this graph. For the spin-singlet the mass shift reads 

AM,.,5^o = ^[-1^(0)1 +2/(2^7 (2^'^^^)^^^^ -^)m] + Oia) = 

= ma\^--^ J- ^)+0(a6). (3.5) 

^8n^ 8n3 16n3(/ + |)^ ^ ' ^ ' 

The evaluation of the singlet-triplet splitting requires some awkward Dirac-algebra, but 
the final result may be brought into a quite transparent form (where one recognizes this 
expression as the well known spin-spin and spin-orbit interaction, adapted to the present 
problem, cf. e.g. pT 



27ra 4 1 f^ — 3|ramP ?>i{r x p)(a*^x dm) 



2 



1 f 2 


3 |tci,7^ I 







m 

Similarly as in the second line of ( |3.6| ) all contributions from graphs 2 are summarized in 
the Breit interaction 

V2 = - —5{r) - " (^p2 ^ r(rp)p\ ^ 

4m3 m? 2m?r \ r"^ ) 

' :LS + ^(^-^-S^]+^S^6ir) (3.7) 
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4 One Loop Corrections 



4.1 One Loop Vacuum Polarization 

In the case of positronium no massless particles are available to one loop order in the 
vacuum polarization, so this effect is only of order of magnitude . In contrast to this, 
QCD contains massless gluons and light quark flavours which may contribute significantly 
to the spectrum. Usually such terms are included in the "running coupling constant". This 
quantity, however, looses its (low order) gauge independence for non-vanishing masses. 
Moreover, only leading logarithms in |g| are contained in that approach. Thus beginning 
with two loops logarithmic terms are "lost" even for massless loops. In accordance with our 
systematic approach we thus evaluate first only the loop in fig. 3. a by standard techniques 
and obtain 



(27r)^r2(g-'f)2 



q'[V -Xnq^ + - -2\n2- (4.1) 



7r2 ^ ^ ^ ■ 3 

7 € 

-e( 2 In 2 - 7 + In 47r) In g ^ + - In^ g ^ + e ■ const + 0(e^)] 

3 2 



with 



1 A- D 

V = --7 + ln47r, e = — — . (4.2) 

We have included the term e ■ f{q) in order to make it applicable in a two loop calculation. 
The graph 3.b contains go terms which can be avoided if we carry out the po integration 
first (cf. (p.28|) and (|2.27|) ). The result can be expanded in powers of f> and a in order 



to show that the effect of go is of 0{a^). With this simplification graph 3.b is exactly 
calculable: 

= |^[iOg2(p _ Ing-^) + 16(7 - 81n2)g2] + 0(e) (4.3) 



Our renormalization prescription consists of a subtraction at the point g = (0, fl), where 
has to be of the order am to avoid large logarithmic contributions from higher orders. This 
seems to be the natural renormalization prescription for bound state problems, because 
also in the BS expectation values (|2.28|), the Bohr momentum am together with po ~ 



0{a m) provides the dominant parts of the integrals [T^, |T^. 



After renormalization, the contribution from the gluonic vacuum polarization (with 



the colour trace already done) reads |25| 



iHg = -«7o®7oi^g (4.4) 
~ 33a\ q' 
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The expectation value of this expression can be obtained by performing the Fourier trans- 
formation into coordinate space, where the integrations can be done analytically (see 
Appendix A). The surprisingly simple result is 

AM, = {Hg) = -ma^^^[^^{n + / + 1) + 7 + In ^] + 0(«^) (4.5) 

IQim^ am 

where is the n-th logarithmic derivation of the gamma function and 7 denotes Euler's 
constant. The closed form of Eq. ( [4.5| ) was not obtained in previous calculations. 
Of course, the contribution to the potential is the standard one: 

= -^(7 + In H (4.6) 

Now we turn to the fermion loops (fig. 3.c). In the literature the lighter quarks are 
usually taken as massless (and 'absorbed' in the number of flavours appearing in the run- 
ning coupling constant) or even ignored |P, ^ , but we will show that they do contribute 



within the order of interest and, furthermore, the explicit dependence on the masses of 
lighter quarks may be important. As pointed out already in the introduction, this is due 



to the fact that the top quark is expected to lie above 170GeV [IS, 18], and therefore 
the bottom and charm quark can neither be taken as relatively massless nor as relatively 
super-heavy as compared to the natural mass scale am. 



The finite part of the self energy in fig. 3.c is a well known quantity |33| for arbitary 



masses mf 



n.^-'^i'Kr-.rn]) (4,7) 
We approximate in the exact solution 



I{q'^,m'^) = J dxx{l — x)\n[x{l — x)q'^ + m"^] 

= ilnm^2-^ + /(p), (4.8) 

fip) = ^+-(i-2p)yrTi;^in^^n±i, (4.9) 

2 

Tfl 

P ■■= ^, (4.10) 



for later convenience /(p) by 
This agrees with the original /(p) better than 1% within the whole integration region. 



/(p)^iln(l + ei). (4.11) 

D p 
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It is easy to transform into coordinate space in order to obtain the potential Vp, 
effectively produced by this fermionic vacuum polarization: 



a 



[Ei( 



5 1 , , 5,, 
- + -ln^ + e3 . 

6 2 mj 



(4.12) 



34 



We note that the energy shift produced by ([4 .121) can be obtained in closed form using 

(4.13) 



r e-^'=Ei{-ax)dx = --[ 
Jo p 



-ax) 



Ei{-{a + P)x)]. 



ln(l + ^)- 

a 



Thus a useful expression for the energy shift induced by fermionic vacuum polarization 
with arbitrary masses mj reads 




'd(3 



2l+l+k 



. 1 , ,^ ^ am 
[--ln(l + /5 ^)] 

P nmjea 



(3^1 



5 
6 



-ln^ + e3 
2 rUj 



with 



(n 



■ k\[{n + iy.f \dp 



p=0 



(4.14) 



(4.15) 



For states up to n = 3 we write this result more explicitly as 

,2 5 , 1 -V 



ma 



And ^ +ln■ 



+ 63 ) 2 



am 



66 + 



(4.16) 



with Ar,i from Tab. 2, 



n 


/ 


am ' 


1 







a 


2 







— + a 


2 


1 




hr' + la^ + a 
- + f - a2 + a 


3 







3 


1 


a' 


- + + \a^ + a 


3 


2 


5 





Tab. 2 
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using the shorthand 



= IH ^ — (4.17) 

am 

„5^ 



Only for mj » am this gives an 0{a ) Uehhng term, modified by off-shell subtraction, 
but for mj — it becomes an 0{a^) contribution, which means that eq. ( 4.16| ) interpolates 
numerically in a range of two orders in a. Therefore ( [4.16| ) or, equivalently, ( [4.12| ) with 
mj 7^ must be definitely taken into account for quarks with mj ^ am at 0{a'^). Thus 
in our present case the finite mass of the bottom quark mf, ^ 0{am) is important. 



4.2 Vertex Corrections 

The one loop corrections to the vertex together with self-energy insertions into the fermion 
lines (fig. 3.d) in the abelian case (positronium), are known to provide corrections only of 
0{a^). The reason for this is the "classical" Ward identity which relates those contribu- 
tions in such a way that the sum of these terms vanishes at |g| — 0. This Ward identity 
happens to continue to hold even in the Coulomb gauge and even in the nonabelian case 
0], but only for the vertex corrections referring to the Coulomb component of the gauge 
field. However, the presence of the gluon splitting graphs 3.e and 3.f produces a contribu- 
tion already to 0{a^) A simple dispersion theoretic argument allows to understand 
this difference: In the abelian case the first graph 3.d in the variable |g| for go = has a 
cut for Re|g| < 2m. Thus corrections in q, e.g. in the electron form factor Fi, for sym- 
metry reasons must be of order g^, because Fi(g^) is regular at 0. This is no longer 
the case with the mass zero intermediate state allowed in 3.e and 3.f. The first - and to 
our knowledge only- computation of the nonabelian vertex corrections in the sense of our 



present approach was performed in ref. |£5|. This work contains certain approximations 



which we want to avoid in order to prepare the ground for future calculations even at 
the level 0{a^). We thus make a systematic expansion and solve the remaining integrals 
analytically which contain contributions of the order of interest. The vertex correction of 
fig. 3.e after performing the colour trace becomes 

SQir'^a'^qi f d'^r 1 / . , rkVi 
g^ J (ivr)^ r^(r — gj"^ r"^ 

x[/^ 7° -/^^ 7l 

7P + 7r — m 7p' — ^yr — m 

^^^'"'^^-(7^..(l,p)7° - 1%M{-I,p'h') (4.18) 



where ( e = ±11 



Vki{e,p) := j-^^— — -^-4* + — 7 — ^ 7^ ^ 4.19 
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After the tq integration it proves useful to proceed with the po integrations (cf. Eq. ( p.28| ) 
), where in contrast to the author of ref. [25|, who approximates aheady at this point, we 
took into account also the pole arising from the denominator of Eq. ( |4.19| ). This results 
in 

. _,/|^z|d_£f(,,,^ ,4.20) 



Fip, r) 

2rEp+,r{r + Ep+,r) - (r + + cu)2] 

(4.21) 

The calculation of the 7 trace to lowest order requires the inclusion of the 7p-terms in the 
wave functions ( p. 21 ) from the projection operators ( p.9|) . After performing this trace, we 
expand in terms of p and ma^ which enables us to combine both terms in eq. ( [4.18| ): 



= J (2^ (f-gT 2r^E. + ^^^^^^ ^'•''^ 

Now Eq. ( [4.22 ) may be evaluated exactly in terms of dilogarithms and the result has a 
cut for |g| < 0, but no pole at |g| = 0. It can be formally expanded for |g| > to 0(|g|): 

The first term in this expansion has been obtained in |]2S|, the second one is the expected 
contribution to O(a^). The BS - expectation value of ( [4.23| ) becomes to 0{a^): 

9ma^ 1 

AMa.e = , ,j^,y 4.24 
From ( |4.23 ) we conclude that the effect of graph 3.e can be summarized in the potential 

Vs.e = ^ (4.25) 

It remains to calculate the vertex corrections with two transverse gluons , graph 3.f. 
At first sight it seems that this graph would give a contribution to order because no 
P7 terms are needed from the wave functions. This, as we will see, is not true because 
the leading (constant) term vanishes as a consequence of renormalization and accordingly 
graph 3.f has been estimated to be of order 0{q.^) in ref. pS. Here we use an approach 
which explicitely provides at least part of the exact result of this contribution. The vertex 
part of the graph 3.f reads: 

?,g^T^ r d^r {qo - 2ro)Y{l{p - r) + m)j'' 
4 J (27r)^ [(p - r)2 - m2](r - g)V 



^ L _ {qi -ri){qk -r-k) _ rjr\ ^ (g^ - ri)r{q- r)rk 
\ * [q — r)2 'P{q — r)^ 
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With the gamma-trace to relative O(a^) and using Feynman parametrization, the effective 
vertex part becomes 

Kfy = /- — TTTT / dy(L-^ -3x^m^L-^)(l + ^ 

'•^ 16 J (27r)^-iio io A;2(g-A;)2^ 

(4.27) 

where 

= {y% " xPoY + k'^ + 2x'pk — 2yqk — x{p^ — m^) — yq"^. (4.28) 

For our estimate it is sufficient to consider the part from the constant term in the second 
bracket in Eq. ( |4.27| ). The k integration can be done easily, leaving us with a finite part 

where 

M = x^m^-w, (4.30) 
w = y{l — y)(f' + x{l — x){p^ — m^) + 2xypq. (4-31) 

This expression cannot be expanded in terms of w because this would yield in the next 
order a spurious linear divergence from the q^ term which would indicate an equally 
spurious O(a^) contribution. Therefore, we expand in terms of {w + y{l — y)q'^) and solve 
the leading part analytically in terms of dilogarithms. Expanding the result in terms of 
(f^, one has 



Vli^ = 4^{2-T[l + ln(^)]^} + other 0(a2) (4.32) 



Since ( [4.32| ) does not contains a term oc the vertex correction due to two transverse 
gluons does not contribute to 0(a''). 

To one loop order also the box graph 3.g occurs in the correction to the BR kernel. 
It possesses an exact counterpart in QED and is known to contribute only to 0{a^) 



35| . Furthermore we have also investigated the two- loop vertex-correction depicted in 
fig. 3.h. Of course this correction is but one of several two loop vertex corrections. The 
renormalization must take into account the whole set of these graphs. Nevertheless, it 
seems that after proper renormalization they yield a contribution to 0{a^). 

5 Other Corrections 

5.1 Two Loop Vacuum Polarization 

As pointed out already in subsection 4.1, the usual renormalization group arguments re- 
lying on massless quarks in the running coupling constant do not consistently include the 
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effect of 'realistic' quark masses in the toponium system, when a systematic BS perturba- 
tion is attempted. However, in the one loop case finite quark masses could yield terms of 
numerical order O(a^), therefore the same can be expected here, leading to corrections of 
0{a^). On the other hand, in a full calculation of effects of 0{a^) two loops with gluons 
cannot be neglected. Although it is enough to consider the two loop vacuum polariza- 
tion for Coulomb gluons only, the computation of all those graphs is beyond the scope of 
our present paper. We just want to indicate how already the graphs 4.a-4.c yield contri- 
butions of this order which are non-leading logarithms. Performing the zero component 
integrations of momenta results in (including a symmetry factor 1/2) 



-iH4,a = -i7o®7o^(n 



(1) 



+ n(2)) 



with 



n(i) 
n(2) 



jD-1 




jD-1 



P 



1 



(27r)^-ip2 



d 



D-1 



P 



q — k — p\ 
[{q — k — p)kY 



(27r)^-iy (27r)^-V2|A;|3|g-'- A;-p1= 



(5.1) 



(5.2) 



(5.3) 



By using dimensional regularization, Feynman-parametrization and usual integration for- 

V{-\ + e)T\l-e) 



mulas we arrive at 

n(i) 



d 



D-1 



P 



1 



r2(l)(47r)i-T(2-2e)^ (27r)^- ^^[(g - 
,(4vr)2^(g2)-2^r2(l-e)r(-l + 2e)r(| 



m 

-t) 



(5.4) 



P(i)(47r)3 



r(|-3e) 



and 

n(2) 



3r(| 



4(4vr)3-2e(| - 2e)(-l + 2e)r(|)r(i)r(| - 3e) 
m\-e)V{2-e){q^Y-^\ 
2(47r)3-2er(|)r(|) 

f r(2e)r(|-6)r(2-6) 

l(-l + 2e)r(|-e)r(f-36) 
2r(2e)r(| -e)r(2 



2\l-2t 



+ 



(5.5) 




4r(2 - e) ^ 4r(3 - e 



r(i-3e) 



From eqs. 



and (|5 



r(|-6)r(|-3e) ■ r(|-3e) 

) the finite, renormalized contribution to the perturbation kernel 



may be extracted as 



A.a 



7o ® 7o 



81a3 ^2 



(5.6) 
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Eq. (|5.6|) differs from eq. (^4.4| ) by a simple factor proportional to a and therefore results 
in the mass shift 

AM4 = + / + 1) + 7 + In (5.7) 

o47r"^n^ am 

The contribution from fig. 4.b is similar. After performing the integrations over the zero 

components we have 



X 



^2 (#)^ i^P)^ , iQP)ipk){qk) 
q — — \ 



This expression can be written in terms of integrals already solved in the course of the 
one loop calculation: 



—r^\V^ + 2V(- 
+2(7 - ln47r - — + 41n2) Ing^ + 21n^ + const + 0(e)]. 



04,,, = - ^ I ^""^ [p^ + 2P(^ - 2 In 2) - 2V In g ' (5.9) 



3 

Eq. (|5.9|) contains overlapping divergences and two graphs like 3. a with one or the other 
vertex replaced by a counterterm have to be added. After that, only an additive infinity 
survives which is subtracted by our usual renormalization prescription. Graph 4.b has a 
net contribution which is proportional to In^: 

The last two loop correction we are considering is the one in fig. 4.c, whose contribution 
to the Coulomb gluon propagator can be written as 



n4.c = -^^^3'j-^j^(<i'^r<i ^ X 



(27r)^ J (27r)^ (g-- kf[q- k - r)^r^k'' 

X + (5-11) 



Eq. ( |5.11|) can be evaluated entirely in terms of gamma functions by a somewhat lengthy 
calculation, but following the same steps as above. The analytic result is 

n • 4 27 2,1-2.,, ,2. r(e)r(2e)r^(2-e)r(|-e)r(i-2e) ^^ 

= -^^2(8^^^) (4^) r(|-2e)r(|-3e)r(l + e) ' ' 

(5.12) 
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Expanding in terms of e and properly renormalizing the result we finally obtain 



n 



4.C 



. 4 q',2Q\ 



(5.13) 



167r4 

In table 3 we collect the results for the self energy parts of fig. 4. a to 4.c, apart from a 



factor — il^g^ 



graph 


n 


4.a 
4.b 


4 /I-^ 

llnC-lln^C 

3 /t^ 2 /I^ 


4.C 



Tab. 3 

The full calculation of the gluon self energy to two loops seems to be very involved in the 
Coulomb gauge. However, already our partial results show the importance of nonleading 
logarithms. 

Two loop calculations of the vacuum polarization are even more difficult if massive 
fermions are included. However, in view of the results from the one loop calculation with 
massive fermions we may expect that nonvanishing masses tend to decrease the importance 
of such terms in practice to something that would be de facto numerically equivalent 
to 0{a^). We try to circumvent these problems, for the time being, by the following 
argument, which also includes the three 'massless' quarks u,d,s. Because of the Ward 
identity for the Coulomb- vertex, it is clear from the theory of renormalization group that 
the same corrections can be obtained by expanding the running coupling constant with 
a two loop (gluons+u,d,s) input for the latter which provides also the first 'nonleading' 
logarithmic contributions. Our present calculation in any case illustrates on the one hand 
the procedure to be followed in a systematic BS perturbation theory. On the other hand, 
we believe that especially the computation methods for the notoriously difficult Coulomb 
gauge may be useful elsewhere. 

The beta function to two loops is renormalization scheme independent for massless 
quarks and its two loop part has been calculated some time ago [p7|] : 



1 

(4^ 

1 

(47r)4 



:ii 



:io2 
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(5.14) 
(5.15) 

(5.16) 



Here is the number of effectively massless flavours and j3{g) is the solution of 

ra dgA 



In 



(5.17) 
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which reads up to two loops 

T ri 1 ^2(/?o + /j,^^) ] 

— o + — In o.^ . ^ -nx • (5.1^ 



/i 2/3o 



Considering this as an equation for g = g{q'^) we 'undo' the renormahzation group im- 
provement by expanding with 'small' g"^ (x a (cf. eq. ( |2.(j| ) ): 



/^9n I 33 — 2nf , (7^ , , 

= a — ^ln^+ 5.19 



Clearly the one loop term agrees with the calculation in subsect. 4.1 in the limit rrij 0. 
That limit, however, is not appropriate here, as stated above (cf. sect. 4.1), because we 
would loose in this way terms of numerical O(ma^). For the computation of the rest we 
need the expectation value of (In^ ^)/^^- This integral can be done analytically (Appendix 
A) and the result is: 

In^ — 2 

(^) = ;^{^ + vl/2(n + / + l) + s„, + [vI/i(n + / + l)+7 + ln^]2}(5.20) 



with 



2(n-/-l)!"g2 (2/+1 + A;)! 



{n + l)\ ^0 k\{n-l-l~ky' 

With eq. (|5.2CI| ) we obtain for the mass shift, induced by the leading logs of the two loop 
vacuum polarization of the Coulomb gluon a contribution: 

^M^ioov = 2 9 27^[— + ^2(n + / + l) + 3n, + (^i(n + Z + l)+7 + ln^)^] + 

1287r"^n^ [ 12 am 

+ 288(^i(n + / + l)+7 + ln-^)l. (5.21) 

am J 

In this expression we have set = 3 as dictated by the number of sufficiently light quarks. 
Whether eq. ( ^.211) really represents the full two loop quark- gluon vacuum polarization, 
numerically consistent with other terms 0{a^), must still be checked in a calculation of 
the Coulomb gluon's self-energy to two loop order in the Coulomb gauge, i.e. going beyond 
the sample calculation here. We, nevertheless, indicate the corresponding potential 

V2100P = -2^y|^ |(33 - 2ri/)'[y + 2(7 + In /ir)2] + 9(102 - yn/)(7 + ln/ir)|(5.22) 
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5.2 QCD 2-Loop Box Graphs 

It would be incorrect to extrapolate from the QED case the absence of corrections to 
0{a^), other than the abelian tree graphs because gluon splitting allows new types of 
graphs. Our first example of a QCD box graph is fig. 5. a. Between the nonrelativistic 
projectors of the wave functions the perturbation kernel from this graph can be written 
effectively as 



dH d'^k {ko + + m) (g° - + - m) 

(27r)4 (27r)4 [(^^ _ ky - m^][{p2 - - m^]k^{q- k^t^iq- t)^ 

1 I (5.23) 



(t - *:)^ \ " ' [i- k)- 

After performing the integrations over the zero components t° and k^ the external mo- 
menta can be set to (m, 0). This is justified a posteriori by the finiteness of the remaining 
terms. The resulting expression will thus only depend on q and m. The leading contribu- 
tion should be 

H5.a = C5.a^, (5.24) 

but a really reliable estimate or even a calculation of the coefficient c^^ is not available yet. 
Supplementing the usual three powers of a from the wave functions for the computation of 
energy levels, we see that the graph 5. a indeed would give a contribution O(a^), expressible 
as yet another correction a — a{l + cs.aCt^) in V which certainly is not included e.g. in 
the running coupling constant. The qualitative result ( |5.24| ) had been noted already in 



|25|, and [Q. It should be noted, however, that also e.g. the graph 5.b may yield a 



contribution of the same structure. A similar graph with crossed Coulomb lines (fig. 5.c) 
is irrelevant because its group theoretic factor vanishes. 

The 'X' graph in fig. 5.d can be checked more easily for possible new contributions. 
As in the calculation of fig. 5. a it can be simplified to give 

^. 6 f d^t d'^k {ko+p\ + m){f + pl-m) 

—iH^d = oiq j - — — r- — — — X 

^ J {27:)^ {2txY [{k + p^Y - rn^Wit + p2f - m^]k\q - kf 

U + Bfni^) . (5.25) 

The integration over k yields a divergence 1/|^ if q and p tend to zero. Contributions 
within the order of interest may only result from possible poles after the t-integration. 
For simplicity we consider the part of eq. ( |5.25| ) from the factor one in the second line: 



dH t° + p° + m 

(27r)4[(t+p2)2-m2]t2(g_t)5 
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rl rl—x 

dx dy 
Jo Jo 



{Any Jo Jo {yq — xp)"^ — yq'^ — x{p^ — m 



2\ 



im /"i , /"i^^ , X 



dx / dy — — - + 0(a) 



(47r)2 Jo "/o x'^rn? + y{l — X — y)q'^ 

In M + 0(a). (5.26) 



(47r)2m m 

Thus the part of graph 5.d, specified above, has a leading term from 

H - 1n 

d2(47r)^|g]m m 

as g — s> and therefore contributes to 0(a^ In a). The second part of graph 5.c gives a 
similar result with a different numerical factor. 



We conclude that - opposite to the QED case box graphs may be important to 
O(a^). They could, in fact, be responsible for changes of the zero order coupling of the 
Coulomb gluon. Substantially different values of that coupling seem to be required in 
phenomenological fits of lighter quarkonia. 

As illustrated by the explicit calculations above, to O(a^) beside abelian QED type 
corrections [^, a host of further non-abelian contributions can be foreseen. 



5.3 QED Correction 

As a rule, the consideration of electromagnetic effects in QCD calculations is not necces- 
sary, but at high energies the strong coupling decreases, and in the case of toponium we 
expect a1 to be of the same order as uqed- We may obtain this contribution by simply 
solving the BR equation for the sum of a QED and a QCD Coulomb exchange results in 
the energy levels 



Po = M° = 2m^l - i^^±^g^ ^ (5.27) 
~ 2m -m— — m— — — + 0(a). 



4n2 2^2 64n'' An? 



where Q is the electric charge of the heavy quark, i.e. 2/3 for toponium. Clearly even the 
'leading' third term can only be separated from the effect of the second one to the extent 
that a(/i) and aqEoilj) can be studied separately with sufficient precision. 

5.4 Weak Corrections 

While also weak interactions usually can be neglected in QCD calculations, this is not 
true in the high energy region, because the weak coupling scales like ^Gpm? ^ becoming 
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comparable to the strong coupling if the fermion mass m is large. Even bound states 



through Higgs exchange are conceivable |38]. Therefore, we have to consider weak correc- 



tions and especially the exchange of a single Higgs or Z particle, assuming for simplicity 
the standard model with minimal Higgs sector. 
The Higgs boson gives rise to the kernel 

- iHHiggs = -iV2GFm^ ^ ^ iv^GFrn"^ (5.28) 

— mjj + mjj 

in an obvious notation. 

Since we do not know the ratio am/niHi which would allow some approximations if 
that ratio is small, we calculate explicitly the level shifts by transforming into coordinate 
space. As in Appendix A, we express the Laguerre polynomials in terms of differentiations 
of the generating function, do the integration and perform the differentiation afterwards 
to obtain 

A , , Gp'm^Oi ^ , am , 

^Mmggs = -m ; , 5.29 

valid for arbitrary levels and Higgs boson masses, with 



-.21+2 n-l-l 



As an illustration some explicit results for the lowest levels are given in tab. 4. 



n 


1 


2 


2 


3 


3 


3 


I 








1 





1 


2 


lnl{0'n) 




a|(2+a2) 




al{Z+&ai+4) 




06 


(l+ai)2 


4(1+02)4 


4(1+02)4 


9(1+03)6 


9(1+03)6 


9(1+03)6 



Tab. 4 



It is evident that Eq. (|]29D will give a contribution of order Gi?m^a^ if the Higgs 



mass is comparable to the mass of the heavy quark and should therefore be taken into 
account in a consistent treatment of heavy quarkonia spectra and the related potential to 
numerical order O(a^). We will thus consider also a corresponding term 



Vmggs = -V2GFm'^^ (5.31) 

correction in our potential. 
Next we consider the contributions of the neutral current, the single Z-exchange and 
Z-annihilation. For the Z-exchange we obtain 

_ q^iq^ 

Hz = -V2Gpml[Yivf-an,)U— ^[Y{vf - an,)],,, (5.32) 

q — niz 
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with 



Vf = -2Qfsm^Q^, (5.33) 
af = Ti, (5.34) 

where is the eigenvalue of the diagonal SU(2) generator for the fermion /. If / is the 
top quark then = 1/2. Because we cannot expect am to be much smaller that mz we 
use the exact expectation value of the Yukawa potential ( ^.30|) . Qq may be dropped in the 
Z-propagator since this provides at most a correction 0(a): 



AMz = -V2G,ml{J^ , '-)= (5.35) 

Gprnia, n 2 -a, ,^ , am , (am)^ ^ , ^ ^ am 
2v27r 3 nm^ Zm'^^n-^ 2 nmz 

S'^ = 5(5' + 1) is total spin of the quark- ant iquark system. Therefore this expression gives 
rise to a singlet triplet splitting within the order of interest. 

Z- annihilation may also yield a sizeable effect. The corresponding energy shift is easily 
evaluated 



2^2 _^,3 



ma 



AMs=o = —TT^—r^io (5.36) 



2V27r n3 

2 2 

AMs=i = ^ AMs=o^lO-^AMs=o (5.37) 

a j 4m^ — m| 

and also produces a singlet-triplet splitting. It should be noted that the two last contribu- 
tions yield singlet triplet splittings which are as important as the usual Breit interaction 
(p.7|). The corresponding contribution to V is 

Vz = v^G^m|4{^[l-i-(5^-3^)(^ + ^)-(5^-^)] + 
■' ^ 2Txr 2aj mzr m^r^ 

^{r) 11 



ml 



■(7- -5 2)} (5.38) 



5.5 Schwinger-Christ-Lee Terms 

As mentioned in section (^, nonlocal interactions have to be added to the Lagrangian in 
Coulomb gauge. We are not aware of any previous attempt to look whether these terms 
give contributions to bound state problems or to some effective potential. 
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By analogy to the second ref. ||26[ we calculate the vi term to 0{g^ 

. 9 



Vl 

KM 



-9 
1 
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47r|p| 



d^rd^r'd^r"A1{r')K,j{r- f')Kjk{r- f")Al{f") (5.39) 
1 



47r|p|' 



This corrects the gluon propagator by 

1 9ig^ 



J d'x J dhdh'dh"-^^Ki,Kju-^^Z,[J]. 



Z[Q]8J;{x^)8Jl{x2) 16 

In momentum space 5G can be calculated by using dimensional regularization to give 

,^'9g'ct 1 



(5.40) 



g5 g2 q2 



+ 



(5.41) 



which means that we have a mass shift with the same structure as the one transverse 
gluon exchange (cf. sect. (^), further suppressed by two more orders in a. 

Since the second term f 2 also represents a correction to the propagator of the transverse 
gluon, it can be estimated by the same method to contribute only in higher orders of a 
as well. We thus find that both terms can be neglected even including terms 0{a^). 



6 Conclusions 

Both, the large mass and the large width of the top quark provide a new field for rigorous 
QCD perturbation theory: In contrast to the lighter quarkonia, (unstable) toponium is a 
weakly bound system, to be treated by Bethe-Salpeter methods in a systematic manner. 
The large width even further reduces the effects of confinement. As shown first by Fadin 



and Khoze |21| V can be included in the (weakly bound) Green function at the threshold 
in a straightforward manner by analytic continuation of the total energy into the complex 
plane. 

In our present paper we first show that the same trick may be also applied to a 
different zero order equation, the BR-equation. On the basis of this result we describe 
how to obtain the proper potential V for such a Green function, rigorously derived from 
QCD in a perturbative sense. Although a fully 4-dimensional formalism is used, which 
especially also allows the inclusion of off-shell effects, our result allows an interpretation 
as a correction to the Coulomb potential: 

Am^ m? 2m?r \ r'^ J 
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5 1 ^ 

(7 + In /ir) + -^^[Ei(-rmjei) - ^ + -ln(-^ + et)] (6.1) 



Svrr 47rr ~^ 6 2 



9 r 47rr •' mi 3 



+x/2G,m,a,^[l - ^ - -3^)(^ + ^)-(5 - ^)] 

It has been obtained by checking contributions up to numerical order 0{a'^) to real energy 
shifts, calculated independently from the weak decay t h + W . The first and the second 
line simply comprise the analogue of the relativistic abehan corrections. The third line 
corresponds to the one-loop vacuum polarization effects for the gluon and from massive 
fiavours. mj = may be used for u,d,s-quarks, but not for c and b. The first term in 
the fourth line results from the nonabelian vertex correction. Among the terms yielding 
0(q;^) in the other entries of the fourth line not all coefficients are known at present : The 
"H"-graphs (5.a+5.b) and the two loop terms involving a heavy fiavour loop (bottom) 
still await an evaluation. On the other hand, we show that certain electroweak effects 
(QED,Higgs,Z) are numerically important to 0{a^). Their respective contributions are 
listed in the last two lines of ( |6.1| ). Comparing (|6.1|) with QCD-potentials proposed pre- 



viously IT^, |TT|, the most important difference is that the coupling as = 3/4a does not 
contain a "running" part. Instead, terms which usually (with zero mass flavours) are in- 
cluded in that running constant have been written explicitely so that orders of the coupling 
constant are not mixed. This we believe to be important for an approach which guarantees 
gauge-independence up to certain order. Of course, all calculations were performed in the 
Coulomb gauge and V should also be used in further calculations only in that gauge. 

Within a rigorous field-theoretic philosophy it would also be incorrect to add, say, a 
linear term to ( |6.1D in order to describe confinement. At best ( |6.1|) could be supplemented 
by a piece oc {G^)r^ which mimics the tail of confinement effects by gluon condensate 

mm- 

In the derivation of our potential we have not only used the level shifts, but also have 
described in much detail new closed forms for such shifts etc. The reason for that has been 
that on the one hand we hope to have given new useful methods to be applicable also for 
nonabelian 0(a^)-effects. On the other hand certain computations of level shifts may be 
useful in conjunction with semi-phenomenological approaches to the lighter quarkonia. 

This work is supported in part by the Austrian Science Foundation (FWF) in project 
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P10063-PHY within the framework of the EEC- Program "Human Capital and Mobihty", 
Network "Physics at High Energy CoUiders" , contract CHRX-CT93-0357 (DC 12 COMA). 



A Expectation Values 

In sect. 4 and 6 we needed the expectation values of logarithmic potentials between 
Schrodinger wave functions. They can be obtained by 

if the expectation value (r^^^) is known analytically. In the latter the representation 

1 Jn-l-l 

,{p) = hm L—-^-—-(l - zy^'-'e"^^ (A.2) 

of the Laguerre polynomials may be used. This allows an easy evaluation of the integra- 
tions and the remaining differentiations can be done with some care afterwards: 

<^-) - ^^^(^' - -) t („ - , ^ - J t ' - 1 

(A.3) 

Using now the Fourier transformations [p9[] 
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F[-^] = -1^-^ (A.4) 



In'^ 1 TT^ 

F[^] = — [-^ + 2(7 + lnHl (A.5) 



one arrives immediately at eq. (|4.5| ) and (|5.2CI|) , respectively. 
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Figure Captions 
Fig.l: BS-equation for bound states 

Fig.2: Tree graphs (broken lines represent Coulomb gluons, curly lines depict transverse 
gluons, wavy lines represent a general gluon and solid lines stand for fermions). 

Fig. 3: One loop graphs and vertex corrections 

Fig.4: Two loop vacuum polarization 

Fig. 5: Two loop box graphs 



30 




31 




32 




33 



4.a 
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